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We show numerically that a rotating, harmonically trapped mixture of two Bose-Einstein- 
condensed superfluids can—contrary to its single-species counterpart—contain a multiply quantized 
vortex in the ground state of the system. This giant vortex can occur without any accompanyiirg 
siirgle-quantum vortices, may either be coreless or have an empty core, and can be realized iir a 
®^Rb-'*^K Bose-Einstein condensate. Our results not only provide a rare example of a stable, soli¬ 
tary multiquantum vortex but also reveal exotic physics stemming from the coexistence of multiple, 
compositionally distinct condensates in one system. 
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I. INTRODUCTION 

According to the conventional paradigm [1, 2], the 
ground state in a rotating superfluid will involve only 
singly quantized vortices (SQVs). Vortices with larger 
quantum numbers are energetically unfavorable and do 
not occur—not even for rapid rotation, which instead 
spawns a triangular Abrikosov lattice of SQVs [3]. Al¬ 
though this is well established [4-6] for a solitary super¬ 
fluid described by a single C-valued order parameter ff, 
vortex physics becomes much more diverse when multi¬ 
ple mutually interacting superfluids are rotated simulta¬ 
neously in the same container. 

Already for the simplest mixture, which consists of two 
superfluid species and is described by two C-valued or¬ 
der parameters and 4'2, a myriad of unusual ground- 
state vortex structures have been found in experimen¬ 
tal and theoretical studies [7]. Experimentally, a ver¬ 
satile platform to study vortices is provided by atomic 
Bose-Einstein condensates (BECs) [8-10], in which two- 
component superfluid mixtures have been realized using 
two different spin states of the same isotope [10-15], two 
different isotopes of the same element [16, 17], or two dis¬ 
tinct elements [18-26]. The unconventional vortex struc¬ 
tures that were detected in these experiments comprise 
coreless SQVs [10] and square vortex lattices [11]. The¬ 
oretical studies, however, have furnished the two-species 
BECs with many more ground-state vortex configura¬ 
tions than the aforementioned two [27-35]: Predicted 
but hitherto unobserved ones include serpentine vortex 
sheets [36], triangular lattices of vortex pairs [35], and, in 
a pseudospin-1/2 representation, giant skyrmions [29, 37] 
and meron pairs [27, 28]. 

One peculiar feature of the two-species mixture, which 
goes against the traditional paradigm, is the appear¬ 
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ance of multiply quantized vortices (MQVs) in the ro¬ 
tating ground state of the harmonically trapped sys¬ 
tem [29, 35, 37]. So far, the MQVs, also known as gi¬ 
ant vortices, have been predicted only in complicated 
states involving a number of accompanying SQVs and 
a large total circulation, thereby requiring rotation fre¬ 
quencies close to the maximum set by the harmonic trap 
frequency. Consequently, the states have eluded experi¬ 
mental observation and, due to the accompanying SQVs, 
might not be suitable for investigating the rarely encoun¬ 
tered ground-state MQV in a controlled fashion. Be¬ 
sides being exotic and interesting in their own right, 
MQVs could also be used to realize bosonic quantum 
Hall states [38], initiate quantum turbulence [39-41], or 
implement a ballistic quantum switch [42]. 

In this article, we make the ground-state MQVs more 
accessible to experiments by showing theoretically that 
an interacting mixture of two dilute superfluids, when 
rotated at moderate speed, exhibits ground states that 
contain a solitary MQV in one of the superfluids. We 
find such states both for mutually attractive mixtures, 
where the MQV has a completely empty core, and for 
mutually repulsive mixtures, where the core is occupied 
by the other, vortex-free superfluid species. These states 
represent a rare instance of a stable, solitary MQV in 
an atomic BEG and, as such, constitute a robust, well- 
isolated, and tunable environment for the experimental 
exploration of MQV physics, in complement to earlier 
observations in mesoscopic superconductors [43-46]. 

All the discovered states share the property that the 
two superfluid species carry unequal numbers of circu¬ 
lation quanta under the same external rotation. This 
requires the two superfluids to be composed of particles 
with sufficiently different masses [47]. For concreteness, 
we will focus on the harmonically trapped two-species 
BEG of ®^Rb and because it has already been re¬ 
alized in several experiments [18-23], it enables a flex¬ 
ible control over its interaction strengths [21], and it 
has a suitable atomic mass ratio of ~2. Although we 
present ground states only for this particular system, the 
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essential features of our results apply generally to mass- 
imbalanced binary mixtures of dilute superfluids. 

II. MODEL 

We assume that the two-species BEC is rotated with 
angular velocity flz. In the zero-temperature mean- 
field regime, the ground-state order parameters dli (as¬ 
signed to ®^Rb) and 'I '2 satisfy the coupled time- 

independent Gross-Pitaevskii equations in the rotating 
reference frame [7]: 

(■^1+5111^11^+5121^3-^1^ - Mi) (g<^) = 0, (1) 

where j G {1, 2}, 

( 2 ) 

and the chemical potentials ensure that J\'i’j\‘^d'^r = 
Nj. Here Nj, rrij, and ojj denote, respectively, the to¬ 
tal number, the mass, and the radial harmonic trapping 
frequency of species j atoms. We only consider quasi- 
two-dimensional configurations pertaining to, e.g., highly 
oblate (prolate) traps with strong (weak) axial confine¬ 
ment and 'kj approximately Gaussian (constant) in the 
axial direction. The intraspecies interaction strengths gjj 
are assumed to be positive, whereas for the interspecies 
parameter gi 2 we also consider negative values. 

We parametrize the interactions by the three dimen¬ 
sionless quantities U = {gn + gi 2 )'m-iNi/h^, 522 / 511 , 
and r = gnj^Jgwgxi > —1- The ground states, i.e., 
the lowest-energy solutions of Eqs. (1), are then uniquely 
specified by these three and the following four other pa¬ 
rameters: 7712 / 7711 , Ai 2 /iVi, wf/wi, and H/wi. Focusing 
on the ®^Rb-^^K BEG, we fix 7772/7771 = 0.471. Equa¬ 
tions (1) are solved numerically using link-variable dis¬ 
cretization [48] and gradient descent. 

III. GROUND-STATE MULTIQUANTUM 
VORTICES 

In order to understand why MQVs emerge in the ro¬ 
tating two-species BEG, we begin with a scenario where 
only r is varied while the other parameters are held con¬ 
stant. Furthermore, for Kj G Z, let (ki,K 2 ) denote a 
sufficiently pointlike phase defect about which arg(il'i) 
winds by ki x 27r and arg ( 4 ^ 2 ) by K 2 x 27r. For interspecies 
repulsion (E > 0), the simplest MQV state appearing as 
the ground state is a (2,0) vortex, whereas for E < 0, the 
simplest one corresponds to (2,1). Below, we investigate 
these two cases separately. 

Consider first the mutually repulsive mixture. Fig¬ 
ures l(a)-l(d) depict ground states at different E > 0 
for a rotating ®^Rb-^^K BEC in which there are two cir¬ 
culation quanta in Rb and none in K. Figure 1(e) shows 
the relevant energy terms as a function of E. When E = 0 


[Fig. 1(a)], the two off-centered (1,0) vortices are sepa¬ 
rated by a distance of ~10 times their core radius. As 
E increases, the condensates move apart, with Rb shift¬ 
ing outward and K inward; this behavior is manifested 
in the trap potential energy, which increases for Rb and 
decreases for K. Consequently, Rb is depleted from the 
region between the two (1,0) vortices, enabling them to 
merge into a (2,0) vortex without the kinetic-energy in¬ 
crease typical of MQV formation; indeed, the kinetic en¬ 
ergy Ti of Rb decreases with E G [0,0.6]. Hence, for 
E > 0.5, we observe an axisymmetric (2,0) vortex, about 
which arg('ki) winds by 2 x 27r. It is a coreless vor¬ 
tex [27-29, 49-52] in the sense that the total atomic den¬ 
sity 77tot = -I- |4'2p does not vanish at the phase 

singularity. We stress that the (2,0) vortex is a unique 
example of a ground-state MQV in a purely harmonic 
trap that occurs as a solitary topological defect without 
any accompanying SQVs. 

The emergence of the ground-state (2,1) vortex for 
r < 0 is illustrated in Fig. 2. For uncoupled conden¬ 
sates [Fig. 2(a)], there are two off-centered (1,0) vortices 
and one central (0,1) vortex. As E approaches —1, the 
two (1,0) vortices move closer to each other, so that 
at E = —0.98, all three phase singularities lie at the 
origin and make up an axisymmetric (2,1) vortex. To 
explain the movement, we note that the kinetic energy 
Ti increases when the two vortices approach each other, 
whereas the interspecies interaction energy £’12 decreases 
due to the increasing overlap f 771772 d^r. Since £12 gains 
in importance when the attraction becomes stronger, it 
eventually begins to dominate over Ti, and thus the (2,1) 
vortex forms [Fig. 2(d)]. 

The ground-state (2,1) vortex constitutes a rare in¬ 
stance of a stable MQV with a genuinely empty, self- 
supporting core. Typically, such vortices are rendered 
unstable against splitting by quasiparticle excitations 
that are highly localized within the core [53-65]. In our 
case, however, the doubly quantized vortex in 'ki is held 
together by the indivisible SQV in '^ 2 - 

The (2,1) vortices are most readily found for relatively 
small values of U. This is because small U implies a 
large size of the vortex cores, which suppresses the kinetic 
energy near the phase singularities and leads to strong 
dependence of £12 on the vortex positions. 

To produce the MQVs of Figs. 1 and 2, it is desirable 
to have control over the parameter E = 512 /^^ 511522 - 
In experiments, gjk may be tuned with Feshbach res¬ 
onances [66], which have been demonstrated for ®^Rb- 
®^Rb [67-69], 41K-41K [70, 71], and s^Rb-^^K [21] in¬ 
teractions. However, ground-state MQVs can also be 
obtained in the ®^Rb-^^K BEC without employing Fesh¬ 
bach resonances. To demonstrate this for an axially uni¬ 
form system, we use the bare s-wave scattering lengths 
Qii/aB = 99 [72], a 22 /aB = 60 [73], and ai 2 /aB = 
163 [18], where ob is the Bohr radius, and accordingly 
set 522/511 = 1-29 and E = 2.27. The remaining parame¬ 
ters are fixed after the ®^Rb-^^K experiment of Ref. [23]. 

Figure 3 shows the resulting ground states at two dif- 
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Figure 1. Emergence of a ground-state two-quantum vortex 
in a rotating, mutually repulsive two-species BEC. 

(a)-(d) Atomic densities |v[ti|^ and |^ 2 |^ and the complex 
phase of the order parameter in the ground state at the 
indicated value of the interspecies interaction strength F = 
512 /^ 511322 - In all four cases, Arg ('Ii 2 ) = const (not shown), 
(e) Interspecies interaction energy _Ei 2 = 512 /|^'i^' 2 p d^r, 
kinetic energies Tj = d^r/ 2 mj, and trap ener¬ 
gies Vj = mjujj A^r/2 as functions of F for the same 

ground states. The parameters in the Gross-Pitaevskii equa¬ 
tions are set to m 2 /mi = 0.471, ( 722/511 = 4, = 10, 

N 2 /N 1 = 1, fl/uji = 0.4, and U = (gii + gi 2 ) miNi/h^ = 
300. The length unit is Ur = ■\/h/miuJi. 


|Ti|2 (s^Rb) ITaP ("'K) Arg(Ti) 



Figure 2. Formation of a ground-state vortex with specieswise 
quantum numbers {k.\,K 2 ) = ( 2 , 1 ) in a rotating, mutually 
attractive ®^Rb-'^^K BEC. (a)--(d) Atomic densities |4iip and 
|'I' 2 |^ and the complex phase of 4(1 in the ground state at the 
indicated value of F = gvij•^g\\g 22 \ Arg [4(2 (r, 0)] ~ (f) for 
all F (not shown). Here 522/511 = ail/i^i = N 2 /N 1 = 1, 
fl/aii = 0.7, and U = 50. 



Figure 3. Ground states of a ®^Rb-'^^K BEC with interaction 
parameters F = 2.27 and 522/511 = 1.29 corresponding to 
the unmodified scattering lengths in a highly prolate trap, 
shown for two different rotation frequencies. In both states, 
Arg ( 4 ( 2 ) = const (not shown). Furthermore, tui/i^i = 2.12, 
N 2 /N 1 = 0.27, and U = 2800 after the experiment of Ref. [23] . 


ferent rotation frequencies, Vl/uji = 0.7 and 0.8. At 
ri/wi = 0.7 (0.8), the Rb species hosts a central 9- 
quantum (12-quantum) giant vortex surrounded by a ring 
of 13 (16) SQVs. In both cases, the K species is vortex- 
free and occupies the core of the central giant vortex. At 
larger 12, the giant vortex becomes surrounded by a tri¬ 
angular lattice of SQVs; similar profiles have been found 
earlier for rapidly rotating single-component BECs in an- 


harmonic trap potentials [74-82]. 


IV. PSEUDOSPIN TEXTURES 

In this section, we analyze our results in the 
pseudospin-1/2 representation [7, 27, 28, 83, 84]. At 
points where ntot 7 ^ 0, we define the local unit-length 
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pseudospin 

® X! (r^k'Hk {r, </>), (3) 

'T-tot 9) -f, 

where cr is a vector of the three Pauli matrices. Now con¬ 
sider a {ki, K2) vortex about which the atomic densities 
are locally axisymmetric. After transforming to shifted 
polar coordinates (r', cj)') with the vortex core at r' = 0, 
we can write for small r', in terms of a spin rotation Z 
and a U (1) gauge transformation acting on a unit-length 
reference spinor y G C^: 

(r', 4>') = (r') ) (4) 

= V^tot ^ Zjk {k^ 4 >') Xk {r '), 

k 

where Ks = K2 — ni and Kg = ki + K2 are integers that 
determine, respectively, the number of 2tt rotations of s 
about the unit vector z and the number of tt windings 
of the U (1) gauge along a contour enclosing the core, 
Z {Kgcj)') = exp{—iKs(j)'(Jz/2), and Cj G M are constants. 

Figure 4 shows s for some of the ground states in 
Figs. 1-3. The (2,0) vortex in Fig. 1(d) is interpreted 
as a doubly quantized skyrmion [29, 37] located at the 
circular interface of the two species, where = |'I' 2 |. 
Since Ks = —2, s rotates by — 47r about z when the inter¬ 
face is traversed azimuthally; additionally, the projection 
z • s changes from — 1 to 1 when the interface is crossed ra¬ 
dially. For the (2,1) vortex in Fig. 2(d), the U (1) gauge 
winds by 37r and the spin s by —27r. However, because 
now z • s vanishes everywhere, this state is not a skyrmion 
but instead corresponds to a singly quantized spin vortex. 
The defect is structurally similar to the so-called cross 
disgyration in the fermionic superfluid ^He-A [85, 86]. 
Finally, the states in Fig. 3 feature giant skyrmions with 
(a) Ks = —9 and (b) Ks = —12. 

V. CONCLUSION 

In this study, we have demonstrated that two-species 
BECs in rotating harmonic traps are able to host ther¬ 
modynamically stable multiquantum (or giant) vortices. 
Such topological entities rarely exist in the ground state 
and have thus been elusive in BECs, whereas they are 
observed and useful in, e.g., mesoscopic superconductiv¬ 
ity [42-46]. In the present case, their stability is not 
induced by elaborate external potentials [74-82, 87-89] 
but is an inherent property of the harmonically confined, 
mass-imbalanced two-species system; The giant vortex 
in the heavier species is stabilized by its coupling to the 
lighter, giant-vortex-free species. 

Experimentally, the presence of the MQV could be ver¬ 
ified, e.g., by measuring the orbital angular momentum 
using surface wave spectroscopy [90-92] or by detecting 
the Kj-dependent concentric density ripples that would 
form in free expansion [93]. Due to its ground-state na¬ 
ture, the MQV is expected to be highly reproducible. 



Figure 4. Pseudospin textures of a (a) two-quantum 
skyrmion [for the state in Fig. 1(d)], (b) single-quantum 
spin vortex [Fig. 2(d)], and (c) nine-quantum giant 
skyrmion [Fig. 3(a)]. The arrows represent the projection 
of the local pseudospin s = onto the 

xy plane. Here cr is a vector of the Pauli matrices and 
ntot = |4 'i|^-I-|4 '2|^. The dashed circle is the species interface, 
where |4/i| = |4/2|. The inset in (c) shows the 2 projection of 
s as a function of the radial coordinate r. 


long lived, and therefore amenable to extensive measure¬ 
ments. 

We also classified the discovered states into spin- 
skyrmion (coreless) and spin-vortex (cored) variants, 
both of which can be realized in a ®^Rb-^^K BEG [18-23]. 
The similarities of these vortices with fractional [94] and 
skyrmionic [95] vortex states in multiband superconduc¬ 
tors, as well as the rich possibilities for the creation and 
tuning of multispecies BECs [18-26], open a wide avenue 
for exploring emergent physics in multicomponent quan¬ 
tum systems consisting of inherently nonidentical com¬ 
ponents. 
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